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Description 

FSELD OF THE INVENTION 

5 [0001 ] The present invention relates generally to object tracking systems, and specifically to non-contact, electromag- 
netic medical systems and methods for tracking the position and orientation of an object. The present invention is 
particularly directed to a novel calibration system and met-c-tr o c r --agn i . ! i i Ir 

can account ft e effec 1 ,\i nonmo* ng metallic objects. 

to BACKGROUND OF THE INVENTION 

[0002] Non -contact methods of determining the pos tic - o si a 

< • that the ect are well known in the art. For example, U.S. patent 5,391 , 1 99, and PCT patent application 
publication WO 96/05768, describe such systems for determining the coordinates of a me Iheter insidi 

the body. These systems typically mck:de one or more colls within the probe, generally adjacent to the distal end thereof, 
connected by wires to signal processing circuitry coupled to the proximal end of the probe. 

[0003] U.S.patent4,71G,708,describesa!ocationaetfc' i ^ ' i , . in, -\ single ax.s solenoid with aferrornagnetic 
core as a radiating coil. There are a plurality of magnetic coii receivers. The position of the solenoid is determined 
assuming that it radiates as a dipole. 
so [0004] PCT patent application publication WO 94/04938, describes a position-funding system using a single sensing 
coii and an array of three, three-coil radiators. The radiator coils are wound on non-ferromagnetic forms. The position 
ofthesensin \ < them 

of the orientation of the sensor coil Is first utilized in order to determine the position of the sensor coll in that order. 
Additionally, the radiator colls ot each array are energized sequentially using a time multiplexing approach. Interestingly. 
?s although this reference discloses that frequency multiplexing can be utilized in order to significantly increase the operating 
speed ofthepositic- system \, h i hu ua \ n lc it complexity. 

It is also important to note that although this reference teaches a single axis sensor position and orientation tracking 
system, it does not address any specific method for calibrating the system. 

[0005] Accordingly, to date, there Is no known system or method that provides for a - - i position sensor 

30 single axis system and method that is capable of being simultaneously driven through frequency multiplexing utilizing 
a novel exact solution technique and a novel calibration method. 

SUMMARY OF THE INVENTION 

35 [0006] Accord 1 c system cc nprising: 

a plurality of field radiators, each field radiator having a plurality of radiator elements, each radiator element generating 
a magnetic field distinct from one another through simultaneous energizing of said radiator elements; 
a signal processor operatively connected to said field radiators; 
«< a test position sensor having three sensor coils: and 

a positioning device; 

wherein said test position sensor Is operatively connected to said signal processor and to said positioning device; 
wherein said positioning device is operatively connected to said signal processor: 
45 wherein said signal processor is pre-programmed with a plurality of test posit thin c rate system defining 

a mapping volur k ; gned with an operating volume produced by the radiators when energised: 

wherein said signal processor Is adapted to guide said test position sensor to each of su d by means of 

said positioning device: and 

wherein sa i \ d ucied to receive a sensing signal from said sensor, indicative of the magnetic field 

50 sensed at said sensor said sensing signal cietining a measured magnetic field at said sensor said signal processor 
also adapted to correlate a know theoretical magnetic field value 1o: each test position: with the actual measured field, 
and to store ( ction of th theoretical magnetic fields. 

[0007] Preferably the positioning device comprises a robot arm. Further preferably, the sensor coils are mutually 
orthogonal. 

55 [0008] According to the invention, there is further provided a calibration method for use with a system capable of 
t ' 3 ( 1 ig a position of a device in accordance with the systems set out above, the method 

comprising the steps of: 
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providing a plurality of field generators, each having a plurality of radiator elements, each radiator element adapted 

to generate a magneiic field distinct from one another througn simultaneous energising of said radiator elements; 

providing a signal processor operatively connected to said field radiators; 

providing a test position sensor having three mutually-orthogonal sensor coils; and 
5 providing a positioning device, wherein said test position sensor is operatively connected to said signal processor 

arid to sa:d positioning device and said positioning device is operatively connected to said signal processor. 

pre- programming said signal processor with a plurality of tesl oosHions within a coordinate system defining a mapping 

volume substantially aligned with an operating volume produced by the radiators when energised; 

bringing the test position sensor to a specified test position' 
10 taking a measurement of the magnetic field of the Z ornpons j t est p > 

stepping the positioning device such that the Y coil of the test position sensor is at the point previously occupied by 

the Z coil; 

a measurement i e r a 3 e i ji p • e ng the Y coil of the test po ersse 

stepping the positioning device such that the X coil of the te position sen iesa ier. intpi iouslyt . c 

• by the 2 coii dunng the first measurement and by the Y coil during the second measurement; 

taking a measurement of the magnetic field of the X component using the X coil of the test position sensor; 
correlating known theoretical magnetic field values with the actual measured fields; and 
storing correlation results for correction of the theoretical magnetic fields. 

so [0009] Preferably, the method further comprises the steps of moving the test position sensor to a next predetermined 
test position; and 

repeating the steps of taking measurements and stepp ig 1 ig device. 

[0010] The present 15 i! I J 

magnetic tleid for tracking a position te device in aooordan e it! the > ferns £ si nt above which method accounts 
25 for the effects of stationary metallic objects that ai p ocati d n and orientation 

medical system is in use. The novel calibration method is used for any medical system capable of generating a magnetic 
field for tracking a position of a medical device. The method comprising the steps of; 

(a) defining a mapping volume within the generated magnetic field; 
30 (b) placing a metallic object within the mapping volume; 

M aligning a sense: it a i _ >m \- ' me mapping Hurt arte measuring me magnetic field at the first point 
with the sensor to establish a first coordinate position (X,, Y h Z,); 

(d) moving the sensor to a next point (X, + dx, Y, + dy, Z. + dz) along one coordinate axis by an added distance 
component (dx, dy, dz) and measuring the magnetic field ax the next point to establish a next coordinate position; 
35 itert metis field at a diafe ' > » f i > 

to establish an interpolated intermediate coordinate position; 

(f) determining the position difference between the interpolated intermediate coordinate position and an actual 
intermediate coordinate position; 

(g) comparing the position difference to an error limit: 

40 (h) setting (X,, Y jt Zj) of the next point as (X, = X, + dx. Y. = Y. + dy. 2, = 2, --- dz) rf the position difference is within 

the error limit and repeating steps (d) - (gj along another coordinate axis; and 

(i) setting the added distance component (dx. dy. dz) by decreasing the value of 1he added distance component if 
the position difference is not within the error limit and repeating steps (d) - (g) along the same coordinate axis. 

45 [0011] The method also includes com g1 nraticn method forthe entire mapping volume in accordance with 
the steps outlined above. Although, the error limit can be any reasonable error range. If is preferable thai the error limit 
be < 1 mm forthe greatest accuracy effects. Additionally, the senso is stepped or move 1 i u 

2 cm to about 3 cm. Moreover, with respect to the stepping of the sensor, the distance moved should remain constant 
oelim-.a ai in h* cal'Drati i Alsc step i ic smpiished by decreasing the value of the added distance 

60 compc ictor of two (X, + dx/ly Y, + d Z, -f dz/:;) 

[0012] A cond embodimen the alibration rr 1c u e vlfh a system capable r metic fieic 

for tracking a position of a device in accordance with the systems set out above, the method accounting for static metallic 
objects comprises the steps of: 

55 (a) defining a mapping volume within the generated magnetic field; 

(b) placing a etaiiic >bject within the mapping volume; 

(o) aligning a sensor at a first point within the mapping volume and measuring the magnetic field at the first point 
with the sensor to establish a first coordinate position (Xj, Yj, Zj); 
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(d) extrapolating the magnetic field of a next point (Xj + dx, Yj + dy, Zj + dz) along one coordinate axis by an added 
distance component (dx, dy, dz); 

(e) calculating the coordinate position at the extrapolated next point based on the pol: in agnetic field to 
establish an extrapc epos lion 

s if determining the position difference between the extrapolated coordinate position and the actual coordinate position 

of the next point; 

(g) comparing the position difference to an error limit: 

(h) setting the; 1 ^mponent(dx, dy,dz)accordingtoaprede - L < i< on difference 
is within the error limst aligning the sensor to a new point within the mapping volume along another coordinate axis 

to j r magnetic fieic point with fhs to establi oint or i \. c 

repeating steps (d) - ;g> along the other coordinate axis: and 

(i) setting the added distance component (dx, dy, dz) by decreasing the ft 6 a d stance component if 

t t u i^j u f e - e i t e - 111 it- „ e| t I ita t[ c it <) steps (d) - (g) 

along the same coordinate axis. 

[0013] The D-eot o mi na sons nt and is preferably approximately 3 cm However, the prede- 

errn 3 i s : c 3 c c J 2 2 L : i tie is il. Addition: s jded distance component can 
be decreased by a factor of two such that the Intermediate point or position is defined as (X, + dx/2 Y, + dy/2 Z, + dz/2). 
[0014] LiC 1 1 1 c 1 t L L t c 1 ih ui 

so according to the vertices of a cube and the entire mapping volume comprises a plurality of cubes. Each cube is defined 
by measurements derived from at least four different vertices Generally, the calibration method is accomplished for a 
mapping volume is approximately 20 cm X 20 cm X 20 cm or (20 cm) 3 . For controlled accuracy In the calibration, the 
sensor is moved by the arm of a robot. 

[0015] These and other objects, features and advantages of the present invention will be more readily apparent from 
25 the detailed description set forth below, taken In conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0016] 

30 

Fig. 1 Is a pempective view of one er dime-it of a - - 1 

Fig, 2 is a schematic illustration of the position sensing components of Fig. 1 ; 

o position and orientation coordinates in 

accordance with a preferred embodiment; 
35 r _ t 1 1 c 1 t l nethod of Fig. 3: 

Fig. 5 is a s I < ~ 1 the present inven 

Fig, 6 is a schematic flow chart showing a method of the calibration system of Fig. 5; 

Fig 7 is a schematic illustration ot an alternative embodiment for a radiator arrangemen having -.adiator coils that 
are not co-located but are mutually orthogonal to each other; 
40 Fig 8 is a schematic illustration ol anofiie alternative embodiment to; a radiate; arrangement having radiator coils 

thai are not co-located and are not mutually orthogonal to each other 

Fig. 9 is a sot d 1 J . 2 J 1 ^ . ■> 1 ■> 1 1 >u : -e . id lie 1 i j '1 t v r 2 y y .tilt 
that are co-located and Put are not mutually orthogonal to each other; 

Fig. 1 0 is a schematic flow chart showing another disclosed method used to find position arid orientation coordinates 
45 with the radiator arrangements of Figs. 7 - 9 in accordance with another preferred embodiment; 

Fig, 11 is a schematic illustration of a calibration cube utilized in a nove c i y xl sr the calibration system 

of Fig. 5; 

Fig 12 is 1 schemata fic»v clml illu jls iin.g orv embodiment of a novel calibration method for accounting for the 
etfects of n n ■> metallic objects utilizing the calibration system 1 Fig. 5: and 
50 Fig 13 is a schematic flow chart illustrating a second embodiment ot a novel calibration method for accounting for 

the effects of nonmovinq metallic objects utilizing ffie calibration system of Fig. 5. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

55 [0017] There is herein disclosed a medical tracking system and method for determining the position and orientation 
of an object, such as a probe or cathetei il 1 1 ntation deten 

method. The system is also utilised with a novel calibration system and method, which comprise the present invention 
[0018] Ot ) and orientation system 20 is shown in Fig. I.The 
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a co >!e 21 wbict enables the use to observ* at functions of a peripN 5 sad as a 

probe or catheter 22 The catheter 22 Is connected to the console 21 at a signal processor 25 (computer) by wire 28. 
The catheter 22 has a single position sensor 26 fixed near the distal end of the catheter 22 along the longitudinal axis 
of thecatheter22nht tion ser r. ratively connected to t e or 25 by wire 28. Console 

s 1 n f> t the signal processc computer) 25 contains signal processing circuits which are typically 

self-contained inside the computer 25. A display 27 and keyboard 29 are operatively connected to the signal processor 
25 wherein keyboard 29 is connected by wire 29a The display 27 permils real-lime display of the position and orientation 
ofthecathet i • 1 "~ i put . ul i f ing in i itstc^irrages 

and Information of the position sensor 26 along with an image or model of particular anatomy such as an organ. One 

to particular use of the system 20 is to map a heart thereby creating a 3D anatomical map of the heart. The position sensor 
6 sedtcc ith a physiological sens as a trade for r pin ion e.g. LAT, 

[0019] Tne signal processor 25 typically receives, amplifies, filters and digitizes signals Tom catheter 22 'Li; 
signals generated by position sensor 26 whereupon these digitized signals are received arid used by the signal processor 
25 to compute the position and o~icnh 0 at the position sensor 26. 

is [0020] The system 20 may also include a reference device 20 « hi _ n^t hown), for 

establishing a frame of reference for the system 20. The reference device 30 is an external patch removably adherabie 
' tf ; exts c ..if'.' 5f a patient artd jperativeiy con d 1o ;h ' x sc 25 by wire 31 . It is important to 
note that the reference device may consist of other alternative forms such as a second catheter or probe with position 
sensor for internal placement within a patient 

[0021 ] Reference is now made to Fig. 2 which schematically illustrates the components of the position sensing system 
20 responsible for determining the position and orientation of the catheter 22 at the position sensor 26, As shown, the 
position sensor 26 is In the form or a coil or hollow cylinder. However, other forms for the position sensor 26, such as 
barrel-shaped, elliptical -shaped (to include asymmetrical shapes! etc:., are also contemplated herein if Is also con- 
templated that the sensor coil 26 can optionally include a flu - canoe: trator as • ;cr« Acidifir nail fie sensor 26 may 

?5 ~ t -it r r r i-rthc ■> -\ - ~, n ->..-.- r ~ ~ ~ r -> a s " aal I effect enso 

[0022] As illustrated and mentioned above, the position sensor 26 is in the form oi a sensing coil and is positioned 
near the distal end of the catheter 22 such that the sensor coil 26 is preferably coaxial or along the longitudinal axis of 
the catheter 22. As defined herein the position sensor 26 can be optionally referred to as a position sensor, location 
e or jsrtiori 1 dc err. - r cc ens >il. coii e term efeis ithete 22 he a 

30 lumen 24 which extends the length of catheter 22. Preferably, sensor coil 26 has a hollow core and is positioned within 
catheter 22 such that the axis of sensor coii 26 lies within lumen 2-1 and along or parallel to the axis of catheter 22. This 
construction provides access from the proximal end of the catheter 22 to its distal end through the hollow core of sensor 
coii 26, permitting the catheter 22 to be used as a delivery device for the delivery of any type of therapeutic through the 
sensor coii 26. Sensor coil 26 and lumen 24 permit the catheter 22 to be used as a delivery device for the delivery of 

35 any type of therapeutic or diagnostic agent or modality or implantable device. For Instance, therapeutics, such as 
pharmaceutical or biological agents, factors, proteins and ceils- tissue repair 0I treatment polymers such as glues and 
adhesives; energy moda es. such as opt : i r ;y uitraso a eguid 

delivery of therapeutic ultrasound, microwave antennas and radio frequency (RF) conductors are just some of the 
examples of therapeutics contemplated for delivery through the catheter 22 due to its unique arrangement of the sensor 

40 coll 26 and lumen 24. Likewise, examples of suitable diagnostics for delivery through the catheter 22. include, but are 
not limited to: contrast agents, dyes, marking fluids or substances. Moreover, implantable devices such as miniature 
It i I mi k H i i <. i a esthetics, etc. are also suitable 

with this arrangement, 

[0023] System 20 further comprises radiator elements or colls 32. 34, 36, 38, 40, 42, 44, 46, and 48 {also known as 
45 generatorcoils). Inoneembcdimenithecoiisarewoundinsetsoithreeorthogor ilani trice forming radiators 

56, 58, and 60, respectively. Preferably, the coils 32, 34, 36, 38, 40, 42, 44, 46, and 48 are each wound around a support 
member such as a spool, in this first embodiment, each radiator 56, 58 and 60 has three coils that are co-located. 
Accordingly, the coils of each radiator are concentric with each other and mutually orthogonal with each other. The 
concentric arrangement is accomplished by having each coil within a particular radiator shaped such that the coils with 
sit their respective support member each have a different diameter. For instance, by way of example with respect to radiator 
56, coil 36 accommodates and receives colls 32 and 34 and coil 34 accommodates and receives coii 32. Thus, coii 36 
(with its support member} has a diameter that is greater than the diameter of coii 34 (with its support member) wherein 
the diameter of ceil Z- 1 h i i 1 il uxrr i I / \ f r h 

arrange 1 all at i - :Js- u .B 56 -.---"0 

ss [0024] in this first embodiment, the radiators 56, 58 and 60 are fixedly mounted to define a plane or location pad 61 . 
The radiators 56, 58 and 60 may be placed in any desired fixed arrangement such as an arrangement in the form of 
i ul alt langie having sides up to 1 meter in length. The radiator coils 32. 34 

36, 38, 40, 42, 44, 46 and 48 are connected by wires 62 to a radiator drive 64. 
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[0025] Generally, radiator drive 64 simultaneously energizes each of radiator colls 32, 34, 36, 38, 40, 42, 44, 46 and 
48 through frequency multiplexing. Each of these colls generates a respective magnetic field (quasi-stationary fieid), 
which induces a voltage in sensing coil 26 which is used as a sensing signal. Accordingly, since each radiator coll 32. 
34, 36, 38, 40, 42, 44, 46 and 48 generates a different magnetic field, the sensing signal that Is generated at the sensing 

5 coil 26 Is comprised of contributions from each of the radiator coils 32, 34, 36, 38, 40, 42, 44, 46 and 48. The sensing 
signal is then fed back to processor 25. The processor 25 then calculates the three position (x, y, z direction) and two 
orientation coordinates {pitch and yaw) (five degrees of freedom or 5 DOF) of the sensor coil 26 from the nine values 
t n i ' c «, ^escribed in detail later below. For instances where the sensing ooi! 26 

has an asymmetrical shape, 6 DOF are calculated to include the roll orientation. 

io [0028] Additionally, there are other alternative radiate * >><.<> r « ' hat a 'e particularly useful with the 
present invention. As shown in Fig. 7, a second radiator a ang mentinc adiat having adi c 

coils 32, 34. 36. 38, 40, 42, 44, 46 and 48 respectively (same coils as mentioned above). The radiator coils 32. 34, 36, 
38, 40, 42, 44, 46 and 48 are mutually orthogonal to each other within their respective radiator 56a, 58a and 60a. 
However these radiator coils are not concentric and are not co-located. Eur rattier the radiators 56a, 58a and 60a 

is comprise a spaced or nonet- e | 5 oup" in which the distance between the colls of each nonco-located triplet 
grouping is no more than three to tour times, the size of the radiator coil diameter. 

[0027] A third radiator arrangement embodiment is shown In Fig. 8. In this embodiment, radiators 56b. 58b and 80b 
have radiator coils 32, 34, 36, 38, 40, 42, 44, 46 and 48 respectively that are not co-located, not concentric and are not 
mutually orthogonal to each other. The orientation of each coil is arbitrary with the only limitation being that one coil is 

;=;; nofp-m f~ - ~r 'he - m ' -nn • / , " - «. K 

[0028] Fig 9 illustrates a fourth embodiment of a radiator arrangement, in this embodiment, radiators 56c. 58c and 
60c comprise a co-located arrangement for the radiator coils 32. 34, 36, 38, 40, 42, 44, 46 and 48 respectively wherein 

5 f j i ji u i i i ntmd mentof "~ig 

the coils of each radiator 56c, 58c and 60c are are not orthogonal to each other. Again, the only limitation to the coil 
25 orientations is thai one coil is not parallel to another coil in a particular radiator arrangement 56c, 58c and 60c, 

Position and Orientation Method 

[0029] Fig. 3 is a schematic flow chart illustrating a method and associated algorithmic components for determining 
30 the coordinates of sensing coil 26, in accordance with a preferred embodiment. The general method steps will be 
described below and the specific steps of the algorithm wiii be described in detail later in this disclose re. Signal processor 
25 (Figs. 1 and 2) determines three positions (X, Y, and 7.) and two orientation (pitch and yaw) coordinates of sensing 
coil 26 by the method described herein. Prior to beginning any medical procedure, the system 20, through the signal 
processor 25, has been pre-programmed with a desired degree of accuracy achieved through calibration. This is a 
35 L 3 >(. u i i i ic r he i r ! I ci nental steps or 

iterations of the algorithm are stopped as soon as the change from a previous step is less than 0.001 cm. The last is 
necessary In order to get 1 mm accuracy). Additionally, for each system 20 that is manufactured, the generators 56, 53 
and 60 through their generator coils 32, 34, 36, 38, 40, 42, 44, 46 ano 4'- . - ' i hhi ted by a novel 

calibration system and method, which comprise the present invention and v h J it detail later In 

40 this disclosure. 

[0030] The system 20, such as the embodiment illustrated in Fig. 1 . Is located in a clinical setting such as a surgical 
suite, and the locator pad 61 is positioned in a desired location. One preferred location is to position the locator pad 61 
near the patient, for instance, beneath a non-metallic operating table (not shown! The system 20 is activated and the 
generator coils 32, 34. 36. 33. 40, 42, 44, 46 and 48 of the radiators 56. 53 and 60 are simultaneously energised wherein 

45 each coi -gnetic fields, each having a distinct frequency. Accordingly, nine separate magnetic 

■quasi-stationary; fields are created Due to the pre-fixed arrangement of the location pad 61. a predictable magnetic 
fieid volume (operating volume) covering approximately a 20 cm X 20 cm X 20 cm or (20 cm} 3 volume (based on the 
radiators 56, 58, and 60 being configured in location pad 61 in a triangular arrangement having sides of 40 cm X 40 cm 
X3 7 a 1 . J ,. v. h^. ■> co w eis •' e are., i. ' > «. i 3 

50 22. it is important to note thai these dimensions are just one of the illustrated examples of size contemplated Smaller 
and larger volume i c ated 

[0031] The catheter 22 is then placed In the patient and brought into the operating volume and the sensor coil 20 
produces sensor signals indicative of the magnitude of the magnetic field at the sensor coil 26. These sensor signals 
are provided to the processor 25 through wire 28 wherein the magnitude of the mag isured field) is 

55 determined and stored in rhe processor 25. 

[0032] As best shown in Fig 3, sit this point the signal processor 25 uses, an i 1 g point (leflecnng both 

position and orientation) and preferably utilizes a dipole approximation (an optional step) for one iteration in order to get 
from this arbitrary point to a starting position for the steepest descent. The arbitrary point Is a pre-programmed starting 
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point in the processor 25. for instance, usually a point chosen at the center of the mapping volume. Then, the program 
in the signal processor 25 moves directly into the calculation with no approximations. However, the convergence can 
> i! ' v ' snsd) by using the dipoie approximation as a first step thus bringing the solution 

closer to the actual position. From there, the full magnetic field with no approximations is utilized. Thus, the arbitrary 

s starting point is used tor the dipoie approximation as well. 

[0033] Once the initial position and orientation estimate is made with the dipoie approximator 72, the processor 25 
calculates the magnetic field at the estimated position and orientation using a magnetic field calculator 74. After this 
calculation, the steepest descent from the calculated magnetic Held of ihe estimated sensor position and orientation 
the measurer; magnetic field of the sensor coil 26 is calculated using a steepest descent calculator 76 which employs 

to the Jacoblan technique. It is important 1o note that the Jacoblan is calculated tor all nine radiator coils 32, 34, 36, 38, 
40, 42, 44, 46 and 48 respectively by the steepest descent calculator 76. The -a r i to 74 and the steepest 

descent calculator 76 use exact theoretical magi 3 i n o da 

[0034] Based on th ilcula - p >n is made, which actually reflects a chanc ! i 

the steepest change in value between tue new estimated position and orientation and trie previous position and orientation 

(the actual position and orientation is not known until the end of the calculation). This result, fax , is added to the 
previous estimate of the sensor position and orientation to arrive at a new estimate of sensor position and orientation. 
Additionally, calibration Information that has been pre-siored in the signal processor 25 in accordance with a novel 
calibration system and method of the present invention is also utilized (greater detail provided below}. Accordingly a 
so pre-determined calibration factor is used to adjust the Jacobian and the fields. 

[0035] At this point, the new estimate of the position and orientation of the sensor coll 26 is compared to the prede- 
termined desired accuracy range. If the new estimate of the position and orientation of the sensor coil 26 is not within 
this range, then ihe steps are repeated beginning with calculating the magnetic field at the new estimated position and 
orientation, 

25 [0038] Additionally, for the altern i'i n n n 

technique is applied by a global converger 77 (Fig. 10) in order to arrive at a new estimate of position and orientation 
that is within the predetermined accuracy range. The specific algorithm Is discussed in detail later below, 

Position and Orientation Algorithm 



[0037] In order to better understand the above-described method, the algorithm utilized by this method will now be 
specifically addressed according to each method step. For illustration purposes, the location and orientation of the sensor 
coil 26 is best described with respect to Figs. 4A - 4C. As shown in Fig. 4A, the center of the sensor coil 26 is positioned 
at point P. The axis of the sensor coil 26 which defines its orientation, is shown as dashed line A in Fig. 4A. 
35 [0038] The position and orientation of the sensor x may be defined as follows. 



wherein Tequais the x, y and z coordinates of the position vector OP as shown in Fig. 4B. The vector n, a measure of 
the sensor orientation, corresponds to the x, yand z coordinates, n x n v , and n z , respectively of A [see Fig. 4C), an 
orientation vector which is a parallel translation to the origin of sensor axis orientation vector A. 
[0039] A partial schematic drawing of the system 20 of the invention is depicted In Fig. 4D in which Q-,, Q 2 and Q 3 
45 are the three radiators 56, 58 and 60 respectively. Each of these radiators 56, 58 and 60 comprises three coils (32, 34, 
36, 38, 40, 42, 44, 46 and 48 respectively). For ease of illustration, radiator Q, comprises radiator coils one, two and 
three {32. 34 and 36): radiator Q ? comprises coils four, five and six (38, 40 and 42); and radiator Q 3 comprises colls 
seven, eight and nine (44, 46 and 48). 

[0040] With this arrangement, through the sensor coil 26, measurements of the fields at the sensor coll 26 due to 
so each radiator coil 32. 34, 36, 38, 40, 42, 44, 46 and 48 Is provided to the proce , c for use w th detailed steps below, 
[0041] Step (a): Determine the initial estimate of sense p^ i i orient u gfl md ienta c 

estimator 70. 

[0042] Fc t^einit 1 ic i and orientation ,e i .it on of the sensor coll 

26 as follows: 



30 




(1) 



40 



55 




(2) 
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[0043] For convenience, A" 0 may be arbitrarily chosen so as to be approximately In the center of the system location 
volume. For example, the three radiators Qi, Q 2 and Q 3 may be viewed as being placed at the vertices of a triangle 
(Fig. 4E). The sides of the triangle may be bisected by medians M-,, M 2 and M 3 . The initial estimate of the three- 
dimensional sensor position / may be chosen to be at the intersection of the medians of the radiators, at a distance, 

20 cm for example, above the plane formed by the radiators Likewise, for convenience. » a may be chosen to be a 
positive unit vector parallel to the z axis. 

[0044] Optional Step (b): Refinement of the estimated sensor position and orientation based on the dlpole approxi- 
mation using the dipole approximator 72. 

[0045] The initial estimate of the sensor position and orientation may be refined using the dipole approximator 72 
based on the dipole approximation, as described below: 

[0046] As shown Irs Fig. 4D, we may define the vectors jj s t R 2 and as the position vectors from the origin to 
the radiator centers for radiators Q,. Q 2 and Q 3 . respectively. S defines the coordinate of the sensor coll 26, For each 
of the radiators, we may define a relationship a, (/= 1 to 3) as follows: f> — j{_ — X 0 „ where X 0 is an Initial position: 

-Mi 



'K\%\ 



(4) 



fl 3=iMi (5) 



[0047] We may define f, for ;' = 1 to 9 as the measured field values at the sensor S of the f elds s sc I 

in h nine i9j. Then ieid I c measured a srd-gto he known techniques 

as out led v. p e t in blic 

[0048] We may also define myfory = 1 to 3 as the sum of the squares of the measured field ; al the s ensor due to the 
colls comprising each of the radiators. Thus, for the system shown in Fig. 2 and Fig. 4D, we have three e 
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[0049] A sensor at points far removed from the radiator coils (distance from point to the radiator > radiator coii radius) 
will experience magnetic fields that are said to be dipolar in character (see for example J. D. Jackson in Classical 
' ,v , . \c-v\ York, 1975, page 178). Under these conditions, as shown 

in U.S. Patent 5,913,820 for the case of three radiators Q 1: Q 2 and Q 3 , each rad ate c * "" r ee concentric and 

orthogonal radiator coils, the distance from the sensor to each of the radiators may be approximated in terms of the 
above-defmed and calculated values of a and m by the following equations: 




vvherein r,, r and ; . are the distances from the sense- to the ;s tit • of radiators Q 1 Q and O, respectively. We may 
use the three distances r P r 2 and r 3 to triangulate to an approximate sensor location / {x,y t z) . Each of the three 
distances r ,. r ? ond ; may be thought of as radii of spheres about each of trie radiators. The iriangulailon process solves 
for the point of intersection of these three spheres, which results in an approximate sensor location as described by 

7{x,y,z) (see Figs. 4B and 4F). 

[0050] Knowing the characteristics of each of the radiator coils (number of windings, coil diameter, etc.) and the current 

passing through these coils, we may calculate fi(l ) , a [9, 3] matrix describing the theoretical field in each of the x, y 

and z directions attributed to each of the nine radiator ceils at the approximate sensor location defined by vector / . 

[0051] Fig 4F shows a single loop radiator coil centered at Q and having a unit orientation vector Q. . The coil is 

centered a distance from the origin of the three-dimensional frame of reference. The vector / connects the origin 
with point M, the approximate sensor location from the triangulation as discurroc x 5 1 • c r : 1 r ^ ,tetec 

to radiator coil / may be found as follows: 

We may define the vector cor such that: cor = / —R t . 

[0052] The vector cor is of n ending t tance I e cc point Q, and point M. 

[0053] We may also define the scalar quantity z as the dot product of O t , a vector describing the radiator coil axis, 
and the vector COT I Z~0 { - COT . 

[0054] The quantity T as shown in Fig. 4F, is equal to the projection of ' cor on O, . 

[0055] We may define the scalar quantity absr as the absolute value of the vector COT \ absr - cor . 
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OS absr'\s equ 1 r and p; 

[0057] We may find B r and B e , vector components of the magnetic flux density at the point Me. 
paraiie! anc el; to the line connecting point M and the x asc din pate 

publication WO 96/05768. 

[0058] We may also def - £ ntities XVBC and xnor as follows: 

xvec = cor - {z • O s ) 



r = |p^|. 



[0059] it follows then from these relationships, as illustrated in Fig. 4F tr 

. . xnor 
sm 8 = 



cos&~ . 

absr 



[0060] We may determine the magnetic field at location M in the frame of reference of the radiator coil ( xi'ec and 
O,) by solving the following rotation matrix: 



- _ (sin 9 cos# Y^'"] _f " 
™ (cos 8 - sin b\b s } ""(<?«,, _ 



[0061] The field f' Q a; :Vl in the frame of reference of the >:, y. ;r coordinate system of Fig 4F dee to a radiator coil 
having a single turn is now given by the following equation: 

40 

{xnor) 



[0082] The calculated magnetic field f at M due to toil ^_ - r v t »n » u it c 3 equation: 



= ^0 ' (effective number of turns of the coil) . 



[0063] 

caused, tor ex I 1 I j lumber of turns of the 

coll may be determined by calibration of the coil, for example, by the method disclosed ZT paten pi 
publication WO 97/42517. 
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006 jcve provides the magnitude of the magnetic field at M due to coi! / in each of the x, yand z 

directions. This calculation is repeated for each of the coiis to provide the data for the matrix H (l) as defined above. 

[0065] We may now define a quantity fo{5£) > the calculated field of a sensor at M corrected for sensor orientation 
n as follows: 



[0066] Making the approximation that the sensor is actually located at the location value given by vector / 5 we may 
then substitute f s the actual measured value of the fields at the sensor for h(x) In the above equation, affording 



[0067] Rearranging this equation affords the following equation which provides a first approximate calculation of the 
sensor orientation vector jj as follows: 



30 [0068] Step (c): Calculate the magnetic field at the estimated sensor position and orientation using magnetic field 
calculator 74. 

[0069] We recall from the Equation above that 
35 A(x) =«(/)•«. 

[0070] Having previously calculated //(/) at approximate sensor location M denoted by vector f, and having 

40 _ 

calculated an approximate sensor orientation n , we may now calculate h\x), the calculated field at position and 
orientation X ~ using the Equation above. The vector is the orientation-corrected magnetic fieid due 

45 to each of the nine radiator coils at a sensor at approximate position and orientation jf = j 7, i? j , and, for the above- 
described system. Is of the form of a [9, 1] matrix. 

[0071] Step (d) C - • die Jacobian) from the calculated magnetic field at the estimated 

sensor position and orienlation to the measured magnetic field at the sensor using the sU i ilculator 76. 

[0072] We may calculate the Jacobian at x - I / , N \ as follows: 



7 /)•*)_ ^(4 
* m, a?,. " 
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[0073] The Jacobian, J lf is the change in the calculated field h(x) at calculated location and orientation x = (/,«] 

for each of the nine radiator coils {/' = 1 to 8), with respect to each of the axc^t riables (/= 1 to 

6), The Jacobian is a 9 by 6 r . - and 6 columns. Ei ;h ■ ement in fr 

jii jch coil / a >sition 3f We then increment ea | ables by a A 

and then sfecv:' ste the f >eid The change In the field due to coil / with respect to each of the six position variables 

p r n ' r ~ k or ~c Tr c r> ^« j 1 "c coils. 

[0074] In practice, :t is quite computationally Intensive to do ail these calculations. Alternatively and preferably, values 
for each of the entries in the Jacobian may be calculated once and stored in a lookup table, individual entries may then 

be obtained by interpolation between lookup table values wherein the increments In x are predetermined In the lookup 
table. 

[0075] We may define the matrix L, a 6 by 9 matrix, as follows: 



[0078] Matrix L is a mean least squares Inversion of the Jacobian. 

[0077] We may compute b , the difference between the actual field at the sensor and the calculated field at the 
25 calculated sensor position and orientation as follows: 

30 

[0078] We recall Equation that: 



h(x) ■£(/)■«. 

[0079] At the actual sensor location, the actual measured field would be equal to the calculated field, which gives rise 
to the following equation: 

[0080] At an approximate sensor location, the field at the approximate location, // (jr) , can be related to the measured 
field at the actual sensor location f as follows. 



7 - £(*__, ) + — - Ax Recall that, 
dx 



*=/-*(*_,). 
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[0082] Thus, we wish to calculate a change ix) in the * i variable X, which 

„ i f I k i t at J ult i l . i c ncuted field at the calculated position and 

1 it -"1 of the sensor. From the parameters and equations defined above, that value ot A xwhich most steeply changes 
the calculated field at the calculated position and orientation in the direction of the measured field Is given as follows: 



[0083] The vector represents the value that is added to each of the six position and orientation coordinates of 
the estimated sensor location and orientation to arrive at a new estimate of the sensor position and orientation for which 
the calculated field is closer to the measured field. 

[0084] Step (e): Calculate new estimate for sensor position and orientation using the new position and orientation 
estimator 78. 



[0085] Having computed a value or Av _ . we may add this value ot ^ to the previous est 
position and orientation to arrive at a new estimate ot sensor position and orientation, as foiiows: 



e of the sensor 



x„ ~ x, + Ax . 



[0086] Additionally, pre -stored calibration information ia calibration factor! according to a plurality ot points within trie 
operating volume is introduced at the magnetic field calculator 74 and the steepest descent/Jacobian calculator 76 In 
order to arrive at the new va!ue/e ia1 lion and orienla 5 the new estimate calculator 78. 

Accordingly, the calit t t i d with the process recited above since it is related to the calculation of the magnetic 

fields at the sensor coil 26 and the Jacobian. This calibration method is described in greater detail below. 
[0087] Step (f): Determine whether the new estimate of position and orientation is within the desired measurement 
accuracy, e g. < 0.1 cm (the accuracy of the system)- However, the incremental steps of the algorithm are stopped as 
soon as the change from a previous step Is less than 0.001 cm as described above which Is necessary sn order to get 
better thanl mm accuracy for the system. 

[0088] One or more criteria may besce.i l » lit t ec. accuracy ot the newly estimated senso; position 
and orientation values. One criterion examines the absolute value of &x the position and orientation correction. If 



PI 



s less than a particular threshold value, it may be a ^ has been calculated 

to the desired accuracy. For example, a value of lAxtelO" 4 era (which results in an accuracy or more than 1 mm) is 

believed to meet the required accuracy criteria for sensors used in most I i - 03tions. 

[0089] The above-described method of determining the sensor position and orientation Involves calculation of the 
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fieid h(^x^ at the estimated sensor position and location x — |/ s Jjj _ We may define the quantity Af as follows: 




[0090] The value oi A, f is another criterion that may be used to determine whether the sensor position and orientation 
have been found to the desired accuracy. For example, a value of Af< 1 0 '• (which results in an accuracy of more than 
1 mm) has been found to be sufficiently accurate for most biomedical applications. 
[0091] 

with the required accuracy. 
[0092] 

position and orientation estimate. The procedure de : bove n steps (c) through (e) is iterated if the est 

a.: • r '"C'ru c < o do not meet one or both of the moc edict acy criteria. Specifically, the newly estimated 
sensor position and orientation values from step (e) are usee; in ctep ic) to recalculate a magnetic -field at the new 
ii j i l j ti n l J luc i„ u u s ep (ch *o refine the steepest 

descent calculation. The refined steepest descent Is used to determine a new estimate for the sensor position and 
c i One fail-safe mech v sm ui red by 'he present invention is to limit the iteration number to a maximum 
number, for Instance, to ten (10) iterations, so that if the algorithm does not converge for some reason, for example, 
due to th n infinite 

loop. The maximum iteration number or iteration number limit is also stored in the signal processor 25. 
[0033] A modification to the ateomrtm described above is employ eel for the alternative radiate i ana smsnts illustrated 
i Figs. 7-9. This modified x n is b id o 1 obai Cc jrgen - 5 tease ned in Numerical Recips I 
052143108 pp383). By utilizing a global convergent technique (as shown in Fig. 10), a more efficient convergence of 

i icntc*i3 r ; achieved Accordingly Ihe 3 
provides for a new correction to the current position AX (as a replacement for the Jacobian step performed by the 
steepest descent calculator 76 of Fig. 3). This global convergence method comprises the following steps: 

First, determine the direction of AX by the formula: 



AX f\&X j. 

Second, find the minimum value of the change In the Field A Fi' Arm in.) along this direction This position is determined 
by: 



AX '~ C* &X (where <XC<1 }. 

Third, update the sensor position according to: 



[0094] It Is important to note that AX correction may not minimize the function if If the estimated sensor position is 
too far from the actual position value. However, the direction of AX is the desired focus at this step. Thus, AFmin. is 

jaio jthe direction of AX. Accordingly, this value wlii be anywhere between 0 and AX. 
[0095] > t tnm doe ^c* cc a erge ( dthm a pred°* <k m mber of iterations; from a certain 

starting position, a different starting point may be selected ^c J^. ofe a(.u t- » i 

e utilized In the global i - ' e is convergence. 

[0098] One helpful technique is to divide the operating volume or working space into subvolumes of 5 cm X 5 cm X 
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5 cm or (5 cm) 3 so that there are a total of sixty-four trials to ensure convergence in any case. Accordingly, this procedure 
may be performed only once. Thus, after finding the first point, the result from the convergent technique is then used 
as an approximation for the algorithm. 

[0097] To summarize then, the disclosed method consists of the following steps: 

5 

Step (a) Estin ate 1 £ ensot position and orientation; 

Optional Step (b): Refine the estimated sensor position and orientation based on the dipole approximation: 
Step (O Z )l i n and orientation 

Step id): Calculate the steepest descent from the calculated field at the estimated sensor position and orientation 
io to the measured field at the sensor; 

Step i«. Calculate a new estimate foi sensor position and orientation which includes utilizing predetermined and 
stored calibration information stored In signal processor 25 in conjunction with steps (c) and (d); 
Step (f): Determine whether the newly estimated position and orientation are within the desired accuracy of the 
measurement; and 

to Step (g): Iterate the calculations (steps (c) - (e)), up to the pre-stoied maximum iteration number, i.e. the iteration 

number limit, at the newly estimated sensor position and orientation to refine the position and orientation estimate 
to the required accuracy (aiso p'e oe i i 3 i 1 . . n ~ 

[0098] Alternatively, Step (d) above is replaced by the global convergent technique described above (Fig. 1 0) (for the 

so radiator arrangements shown in Figs. 7-9). 

[0099] In practice, the disclosed method is applied to the sequential calculation of multiple number of sensors positions 
and orientations at a plurality of points In space. For the calculation of successive points that are closely spaced :n time, 
it may be assumed that the sensor has not moved significanllv icel s previously determined values Accorc i 
good value for the initial estimate of position and orientation for the n ;h position and orientation will be the values 

25 determined at the n-1 st position and orientation. 

Calibration System 

[01 00] A calibration system 90 according to the present invention is operably connected to the position and orientation 

30 j " 2 1 J li il i ' 1 f i 3 i ir J pic i 

in the signal processor 25 which is used in performing an accurate determination of the sensor coil 26 (Fig. 2) position 
and orientation when introduced into the operating volume. As shown in Fig 5 -he c 3 n 90 includes primary 

components of the position and orientation system 20 including the signal processor 25. the: radiator drive 64, the location 
pad 61 and radiator or generator coils 32, 34, 36, 38, 40, 42, 44, 46 and 48 wound around spools forming radiators or 

35 generators 56, 58 and 60 respectively along with the respective wires 62, 

[0101] The calibration system 90 further includes a test position sensor 100 having three sensor coils 102, 104 and 
106 which i \ i 1" < in t in 1 -i c pi J T pit ci -i 

between the colls in the sensor 100 is 1-0.02 cm. The test position sensor 100 is operably connected to the signal 
processor 25 by wire 28 for processing sensor signals provided by the test position sensor 100. Additionally, the test 

40 position sense; 1u0 is ope: abk connected to a positioning de-iee such as a robot aim i 10 by .,v-e i id The robot arm 
1 10 is operably connected to the signal processor 25 by wire 114. The processor 25 has predetermined test positions 
according to the three coordinate system, e.g X v'__xcinit I i „ led and fifty 

(1 50) test positions are predetermined and stored in the processor 25. These test positions closely align with the operating 
volume produced by the generators 56, 58 and 60 when enerc c 3 i e processor 5 5 programmed with each 

45 of these test positions, the robot arm 110 can be accurately guided and positioned at each test position. This method 
will become much clearer when described in greater detail below. 

Calibration Method including Algorithm 

50 [0102] In general, in the calibration method according to the present invent of the location 

pad 61 (Fig. 5) are adapted to the measured field 1 a specific position and orientation system 20. Accordingly, each 
system 20 is customized with its own calibration data which is used during the location calculations when determining 
the position and orientation of a sensor coi! 26 (Figs. 1 and 2} as mentioned above and depicted In the schematic flow 
chart of Fig. 3. Pending US patent application no. 09/180,244, entitled Rac iblished as US 

55 6,335,617, also addresses techniques for calibrating magnetic field generators or radiator coils. 

[0103] Fig. 6 shows a schematic flow chart tor the calibration method according to the present Invention. Accordingly, 
with this calibration method th i \ r ~r -i i 

Although ire 1 oris are preferably utilized, any number of test positions may be utilized 
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i 1 f\s mentioned 

abo th I volume, for instance, 20 cm X 20 on X 20 

cm or (20 cm) 3 . As defined herein, the terms "mapping volume", "sampling volume" and "mapping area" are analogous 
terms for "operating volume" and are used interchangeably throughout this disclosure. 

5 [0104] First the radiator coils of radiators or generators 56, 58 and 60 of the position system 20 are activated through 
simultaneous energszmo using ilrequ nc, multiplexing technique in orde* to create nine distinct • \ ,n ;t; fields (quaa- 
stationary magnetic f-l j - I !rntneoperafingvolume.Thentherobotarm110lsguidedandacci.irateiypositioned 

h re c ' < : Y a: omponents of the magnetic lield are r 1 

position sensor 100. After each measurement, the processor 25 checks the sampling of the j - 1 ens in order to 

io determine if a requisite sampling volume has been achieved. The sampling volume corresponds to the operating volume 
of the field radiators 56, 58 and 60 Once the desired sampling volume has been completed, the processor 25 correlates 
the known theoretical magnetic fields at each test position (X,. < 7 to th actual measured field at each test position 
(X., Y,. _ This correlation :s a mathematical transformation that maps the theoretical magnetic fields at the acquired 
in * \ i ui \ 7 ~~ in i in I i d 

to and the results are stored In the processor 25 foi each pc lull > x Z ■ ano n- - p ^ mJ orientation 

sensing procedure such as that previously described, for example, with the system 20 depicted in Fig. 1 . The calibration 
data from these stored calibration files is used to correct the theoretical magnetic fields. 

[0105] One example of the calibration method - ip r * ation system 90 (Fig. 5) is outlined below. As shown, 

using the calibration system 90, the X. Y. Z components of the magnetic field is measured using the three-axis test 

so position sensor 100 comprising three orthogonal sensor coils 102, 104 and 106 which provide sensor signal information 
to the processor 25 of the po< tr - Mhe mapping between the location pad fields and the 

theoretical fields one needs to know the exact coordinates ol the measurement with respect to each of trie coiis 32, 34, 
36, 38, 40, 42, 44, 46 and 48 respectively of the iocation pad 61. To this end, the robot arm 1 10 Is aligned with the 
location ps i - ; i i \ e 3 e predefined e c ion ' Z le acceptable 

25 sequence is as follows: 

a) Bring the robot arm 1 10 to a so- r,/ f ~ t m \ * X is a test position for measuring relative to the 
Z coil 106 of the test position sensor 100. 

b) Take a measurement of the magnetic field of the Z component using the Z coil 1 06 of the test position sensor 1 00. 
30 c) Step the robot arm 1 10 such that the test position sensoMOO is moved 1 cm forward (this places Y coil 104 at 

the same point previously occupied by the Z coil 106)= and measure the magnetic lield of the Y component using 
Y coil 104 of the test position sensor 100 . 

d) Step the robot arm 110 such that the test position sensor 100 is again moved 1 cm forward (placing the X coii 
102 at the point previously occupied by the Z coil 106 during the first measurement and by the Y coi! 104 during 
35 the second measurement) and rr e the 1 ^ N sr r. n inj / c pi, 

sensor 100 . 

[0108] As mentioned ■> -> - t ^ i e „ pei f r <=o r- j ^nphnq o ^n e of 20 cm X 20 cm X 20 

cm around the center and above the iocation pad 61 which corresponds to the operating volume when the system 20 
40 (Figs. 1 and 2) is used for medical purposes. Additionally, the steps between two successive points/positions are 2-3 
cm in each direction though the steps need not be constant through ail of the sampling volume. Typically, there are one 
hundred and fifty (150) points sampled ' v i .'if m w att e end of the ca static nethod e 1 

{B'j I j = 1 ... 9,i = 1 ... * 150 } 

[0107] Where the i Index is the number of points ^ jc< r 3, s the frequency idex to r r jene 1 1 : il a tis ;1 
frequencies. Thus, the index / refers to the measured fields for all of the test positions, e.g. 1 50. 
so [0108] The mat en of space there exists 

a rotation transformation between the measured fields and the theoretical fields. Thus, we can rely on the equation: 



OB n = B M (6) 

whereby O is a 3x3 rotation matrix that is given by the equation: 
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Q-B^B^ (?) 

5 [0109] It is important to note that both B m and B 1h are [3x150] matrices and the mat'- e al trans ormation maps 
B m and B. h as close as possible in the least mean square sense. Also, note that in equation (7), we assumed that the 
same transformation is applied to the entire space volume. Irs fact, we may have a separate rotation matrix to each 
surovolume o1 our mapping space, For instance, one can o':v:de the mapping volume mio sub-cubes having a sub-volume 
cm X 5 cm X 5 cm. e.g each sut \ of 

io the sub-cubes. This rotation matrix is then stored in the signal processor 25 with the caiifc inf i '<■ on regarding 
its position in space. This procedure is conducted for each radiator coil 32, 34, 36, 38, 40, 42, 44, 46 and 48. Typically, 
one can derive approximately twenty (20) - thirty (30; 3x3 matrices for each coil 32, 34, 36, 38, 40, 42, 44, 46 and 48. 

[0110] Accordingly, if one were to begin at a starting point X, then the theoretical field 1 poir B th (if) can be 
jicuiated by the algorithm of present invention Then. 1 new" theoretics.: field Is given b; thi c xpression: 



0(x)B(x) 

[0111] Where o(x) is the relevant transformation at the point x. And for the calculated theoretical Jacobean to be J(x), 
then the "new" calibrated Jacobean is given by the equation: 



25 0(x)J(x) 

[01 12] A;d r:, x wci B {magnetic field) and J (Jacobean) are used in the algorithm of the present 

invention in the same manner that noncailbrated versions are used. However , with the added calibration information 
so ;the calibrated B). the system 20 has greater position accuracy which is particularly useful in medical procedures 

Calibration Method with Static Metal Compensation 

[0113] T he Dissent invention also includes a novel calibration method tor the radiators 36. 58 and 60 of the position 
35 and orientation system 20 (Fig. 1) capable of compensating for the effects of disturbances caused by nonmoving or 
static metallic objects placed within the mapping or operating volume. The calibration method :s used in conjunction 
with the ca c i 3 including th s a ~ 1 f" n s a - 1 

method is useful for various arrangement of the radiators 56. 58 and 60 to include the radiator arrangement embodiments 
shown in Figs. 5, 7, 8 and 9. 

40 [0114] Accordingly, a metaiiic object (not shown), such as the C-arm of a fluoroscopic device, is placed within the 
intended mapping volume generated by the radiators 56, 58 and 60. The purpose of this step is to establish calibration 
data for storage in the s jnal pr -> i t in u t u t f r i 1 1 r ct 

onthemagni f - e rated by the radiator colls 32, 34, 36, 38, 40, 42, 44, 46 and 

48 of the radiators 56, 58 and 60 respectively. 

45 [01 15] As mentioned above, a preferred mapping volume is the magnetic field inside an area approximately 20 cm X 
en X Ocn or i cm Fhs sal it tion l hod a t { \ , - t 'i" *~ - .a " - " in ~tk 
field produced by any rnetal object located within the mopping volume. Once the magnetic field is mapped according to 
the method of the present invention described in greater detail helow tor instance, using the single- axis sensor algorithm 
associated wit!] the system 20 (Fig. 1). one can use the system 20, even In the presence of one or more metallic objects, 

so as if it was a com e etic field situation, i.e. a situation without the presence of metallic object 

interference. 

[0118] One embod r t od according to the pres best liiustra 

of Rg 1 2. First, an intended mapping volume is def - static metallic obje ; d within the 

intended mapping volume. Again U ould be those cb--. le field of the 

55 patient and within the mapping volume. One example of a typical nonmoving (nonmoving meaning only relative to the 
location p. . a ; u' ^' ' .j \. -p. J; :e -*ccc c igl; the C-arm is placed within the mapping 

volume. 

[01 17] Turning now to Figs. 5 and 1 1 , the robot arm 1 10 moves the position sensor 100 to a starting point to begin a 
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mapping ar J 3 a mapping calibration cube 200a. For e<ar ng point is a 

position or point beginning at the vertex of a cube having a volume of (3 cm) 3 . For example, an appropriate first point 
or starting position is the; point 210 of the cube 200a shown in Fig. 11. 

[01 1 8] The sensor 1 00 :s aligned at the first point orcoordinate position 210 represented according to three dimensional 
5 coordinate (X,. Y, Z,) within the mapping calibration cube 200a (within the mapping volume) and the magnetic field of 
the radiator coils 32, 34, 36, 38, 40, 42, 44, 46 ,48 is read and measured at the first pc sh the first coordinate 

position 210 with the sensor 100 through the signal processor 25. The sensor 100 is then stepped or moved to a next 
or second point 220 (Xj + dx, Y; + dy, Z, + dz) defined by the first position and an added distance component (dx, dy, 
dz) by the robot arm 1 10. The magnetic field at the second point 220 is then read and measured through the signal 
10 processor 25 to establish a next or second coordinate position. 

[01 19] The signal processor is preprogrammed with the calibration positions of the calibration cube 200a which are 
1 >3 > jf a (3 cm) 3 volume cube. The robot arm 110 steps the sensor 100 In 

- ^fcn 2 3 e 1 e e a he stepped d ce s 3 1). The exact 
step distance is predetermined and remains constant throughout t^ mappir c th the only ex- 

;i ; cepfion being vvhere accounting for metallic disturbance effects at a particular position. This modification to step distance 
1 ceiow 

[01 20] For example, for the procedure schematically illustrated in Fig. 1 1 , the step distance is 3 cm (along one of the 
coordinate axes X. Y or Z; tor each vertex point or position 210. 220. 230. 2-10 and 2o0 Additionally, the robot arm 1 '0 
steps the sensor 1 00 according to the direction of the arrows along the edges of the cube 200a (along each coordinate 
so axis), 

[0121] Once the magnetic field is measured anc ositi J trained 0 poin 0. e.g , x emit 

position (Xj, Y|. Zj) and at the next or second point 220, e.g. next or second coordinate position (Xj = X, + dx, Yj = Y, + 
dy, Zj ~ Z, + dz), reflecting the added distance component dx, dy, tiz where dx = 3 cm (also referred to as the step or 
teppeJoi tanoe) t -ir 1 ne - 1 n intermediat 2 e t= p n the first or^„ tinn 210 esec positic 

25 220 is interpolated and the position coordinate of the Intermediate point 220a is calculated by the signal processor 25. 
it is important to note that the intermediate point 220a will lie along the coordinate axis of the stepped distance, for 
Instance, within a 3 cm distance along the X coordinate axis between vertice points 210 and 220 
[0122] After interpolating 11 3 r 1 1 J " " il it (.interme- 

diate point (calculated intermediate position using the position and orientation algorithm}, the signal processor 25 takes 

3D the difference between the calculated intermediate position and the actual intermediate position. The intermediate 
position difference (e) is then compared to an error limit also prestored in the signal processor 25. Although the error 
limit may be any value, an error limit <1 mm has been found to be reasonable and an acceptable error limit. 
[01 23] The next or second point (X.. Y, Z ; ) Is then set as (X, = X, + dx, Y t ~ Y; + dy, Z, - Z, + dz) by the signal processor 
25 if the position difference (e) is within the error limit (<1 mm) and the robot arm 1 10 steps the sensor 100 to another 

35 point 230, e.g. a third or another next point (third coordinate position along another coordinate axis, e.g. Y axis (third 
vertex of cube 200a along step distance dy = 3 cm). 

[01 ?d] lire -1 r ( fa nt eg i s not < 1 mm, the signal processor 

25 decreases the value of the added distance component, e.g. sets dx, dy, dz as dx=dx/2, dy=dy/2, dz=dz/2 . For 
example, the added distance component dx Is dec 1 i 1 n (3 cm : 2 err snd whereby dy dz 0) and the robot 

«< arm 1 10 repeats the step to a new second point 220b along the same coordinate axis, e.g. X axis. Accordingly, the 
sensor 100 Is stepped an added distance component of dx=1 .5 cm to the new second point/vertice 220b whereby the 
magneto "led s ead at pc rt 22Di> u ■■d 1. . 3 i r 1 1 t 1 <. c 

is then interpolated for a new intermediate point 220c, e.g. a point/position between the new second position 220b and 
the first position 210 and the position of the intermediate point 220c is calculated (new intermediate position) using the 

45 position and orientation >•>•■) And, just as before, the new inte me ec position is compared to the actual 

position of the intermediate point 220c to determine if this difference (e) is within the error limit (< 1 mm). If the position 
e ns> »« t n • - - 2 nit the robot arm steps the sensor 100 to yet another point vertice 230a (position 
erik 3 5 cube 20r ilong anoth sordln ix g. Y a> 1 added distance component dy=3 cm and 
repeats the steps outlined above for each vertice c 1 cube 200a. 

50 [0125] As mentioned above, the added distance component of the inte m c f asing the value 

of the a del c ioi within the error limit "She 

the decreased adjustment to the added distance component 

[01 26] This calibration mapping process accounting for static metal compensation is continued according to a series 
or c uralit\ - .-on cubes beginning with the mapping of a second 0 3b until the 

55 complete mapping volume of 20 cm X 20 cm X 20 cm or (20 cm) 3 is completely mapped and calibrated accounting for 
existing met. - aimn the mapping volume, it Is Important to note that each of the at on cubes 

200a, 200b, etc. may very well be cubes having sides thai are not of equal lengths due to the effects of a metallic object 
encountered during the mapping of a particular cube. 
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[01 27] A second embodiment of the calibration method according to the present invention utilizes exlrapolaiion of the 
magnetic field at the next point. As shown in Fig, 13, in particular, first, the meta 1 t is f :ei within the intended 
mapping volume. Again, the robot arm 11 0 positions the sensor 100 at the first point 210 (starting or first position X l; Y., 
Z) and the magnetic field at this first point is read and measured by the signal proce^ determ ne the first 

5 3 r ition. r i r 2 - 31 he magne c fieid at the next or second point 220, e.g. 

next or second cocic ex vertex of mappin ube 20 » c i the exf ipo ed point - ed a X. = 

Xj + dx, Yj = Y, + dy, Z; = Z, + dz) which includes the added distance component dx, dy, dz as appropriate, in this case, 
the added distance component (dx) is a 3 cm distance along the X coordinate axis. 

[0128] On tn rsagnetsc fista for th< cone jint th« " v t oint extrapolated, the location coordinate of this 
io point 220 (calculated second position) is calculated using the position and c v n 1 w >mpared to the 

actua ;os:tioacoe:di < iertodetern ne he pos on difference tdiffe nee 

hi 1 1 i T ' i nsor 1001c 

e' point 230 (next <- 1 ~ along another coordinate ax:s. e.g. Y ax:s and - - steps above are continued along the Y 
coordinate axis, etc. if the ? thin th r lln r i i is exti 

is for an Intermediate point 220b, e.g. intermediate position coordinate (X, + dx/2, Y j + dy/2, Z. + dz/2) between the second 
point 220, e.g. second coordinate position and the first point 210, e.g. first coordinate position (still along the same 

o x' hi At 1 J i eei l ! l l 1 1 n, 

t i ease the sd J c corn n if A rnagn ?id measurement is then taken at 
the Intermediate point 22Gb with the sensor 100 and signal processor 25 in order to determine the intermediate position 
c po a orientation he intermediate posi hich is actua y a new 

second or next point), the rest of the steps arc Ml i 11c i h |. d 1 ff r r 

(e) of the intermediate point 220b (new second point) and determining whether e < 1 mm. 

[0129] The mapping calibration procedure is continued in accordance with the method steps shown in Fig. 13 thereby 
creating r i me is completed (approximately 20 cm X 20 cm X 

25 20 cm or (20 cm)3 volume), 

[0130] It is clearly contemplated by the present invention that the steps of the methods and algorithms outlined 
throughout this disclosure may be sequenced differently than the sequences described above which served for example 
purposes only, it is clearly contemplated that the sequencing of these steps may be alternated and, in some instances, 
optionally deleted and still satisfy the novel requirements of the present invention, 

30 [0131] It will be appreciated that the preferred embodiments described above are cited by way of example and the 
full scope of the invention is limited only by the claims, which follow. 



Claims 

35 

1. A calibration system (30) comprising: 

a plurality of field radiators (56 58 60) each field radiato: having a plurality of radiator eie cents (32-48; each 
radiator element generating a magnetic field distinct from one another through simultaneous energizing of said 
40 radiator elements (32-48); 

r- -v -,| s ~ m r r 

a test position sensor (100) having three sensor coils (102. 104. 106): and 
a positioning device (1 10): 

ition sensor (1 00) is operativefy connected to said sigr 21 ana to said positioning 

device (110); 

wherein sad positioning devsce (110) is opeiatlveiy connected to said signal processor (25); 
wherein said signal processor (25) is pre- programmed with a plurality of test positions within a coordinate system 
defining a mapping volume ub ^ j 1 i c i t i J ji 5 8 6 

so when energised; 

i io 6 A l 2t t ms 

by means of said positioning device (110); and 

wherein said signal processor (25; is adapted to receive a sensing signal from said sensor (100). indicative of the 
magnetic field sensed at said sensor, said sensing signal defining a meas.m •> - r < -n w ; . 
55 signal processor also adapted to correlate a known theoretical magnetic field value for each test position with the 

actual measured field, and to store correlation results for correction of the theoretical magnetic fields. 

2. The syster i s tioning device comprises a robot arm (1 1 0). 
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! he system according to Claim 1 or Claim 2, wherein said sensor coils (102, 104, 106) are mutually orthogonal. 

A calibration method for use with a system (90) in accordance with claim 1 capable of generating a magnetic field 
for tracking a position of a device, the method comprising the steps of: 

providing a plurality of field generators (56, 58. 60). each having a piura 1 32 18 eac 

I sen 5 iagn i-f dddistinctfrom one anotberthrough simultaneous energising 
of said radiator elements; 

provid ->g a s gna p see sc pe co nec 1 >& 60); 

providing a test position sensor (100) having three mutually-orthogonal sensor coils (102, 104, 106); and 
providing a positioning device (110). wherein said test position sensor i 100) is operative;;/ connected to said 
signal processor (25) and to said positioning device (110) and said positioning device (1 10! is operatively 
connected to said signal processor (25); 

pre-programming said ria i plurality of test j t dinate system defining a 

mapping volume substantially aligned with an operating volume produced by the radiators (56, 58, 60) when ener- 
gised; 

bringing the test position sensor i 1 00) to a specified test position (X 0 , Y 0 , Z 0 ); 

taking a measurement of the magnetic field of the Z component using the Z coil (106) of the test position sensor (100); 
stepping the positioning device (110) such that the Y coil (104) of the test position sensor (100) is at the point 
previously occupied by the Z coil (106); 

taking a measurement of the magnetic field of the Y component using the Y col! ( 1 04) of the test position sensor (100); 
stepping the positioning device (1 10) such that the X coil (102) of the test position sensor (100) is at the same point 
previou v occupied by the Z co i j the f - irid b e Y coil (104) during the econd 

measurement; 

taking a measurement of the magnetic field of the X .component using the X coll (1 02) of the test position sensor (1 00); 
correlating known theoretical magnetic field values with the actual measured fields; and 
i r kn I u ri< < h il r i ii 

The method of claim 4, further comprising the steps of moving the test position sensor (1 00) to a next predetermined 
test position; and 

repeating the steps of taking measurements and stepping the positioning device (110), 

A calibration method for use with a system (90) in accordance with claim 1 capable of generating a magnetic field 
for tracking a position of a device, the method accounting for static metallic objects and comprising the steps of: 

(a) defining a mapping volume (200a) within the generated magnetic field; 

(b) placing a metallic object within the mapping \ olume. 

(c) aligning a sensor (100) at a first point (210! within the mapping volume and measuring the magnetic field at 
the first point with the sensor to establish a first coordinate position (X., Y., Z.): 

id) moving the sensoi to a next point (220: X, + dx. Y, + dy, Z, + dz) along one coordinate axis (210, 220, 230, 
240, 250) by an added distance component (dx. dy. dz) and measuring the magnetic lield at the next point to 
establish a next coordinate position: 

(e) interpolating the magnet c o mm t „v » " „v > f p j < t^nd the next 

ooid e w i - sh an - c 1 - tetmediate oordinate position; 

i t i 1 -m <il i t i -i i 1 r 5i i Oai and an 

actual Inten e^ e position, 

(g) comparing the position difference to an error limit; 

(h) setting (X„ Y„ Z) of the next point as (X; = X; + dx, Y = Y, 4 dy. Z. = Z, + dz) If the position difference is within 
the error lin c teps (d) - (g along ai f axis (210 220 230 240 250); and 

^ri^iJ u 1 ciit 

if the position difference is not within the error limit and repeating steps (d) ■ (g) along the same coordinate ax is 

A calibration method for use with a system (90) in accordance with claim 1 capable of generating a magnetic field 
for tracking a position of a device, the method accounting for sta me": 

(a) defining a mapping volume (200a) within the generated magnetic field; 

(b) placing a metallic object within the mapping volume; 
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(c) aligning a sensor (100) at a first point (210) within the mapping volume and measuring the magnetic field at 
the first point with the sensor to establish a first coordinate position (X,, Z,); 

(d) extrapolating the magnetic field of a next point (220; X, + dx, Y, + dy, Z, + dz) along one coordinate axis 
(210, 220, 230, 240, 250) by an added distance component (dx, dy, dz); 

5 (e) calculating the coordinate position at the extrapolated next point (220) based on the extrapolated magnetic 

field to establish an extrapolated coordinate position; 

(f) determining the position difference between the extrapolated coordinate position (220) and the actual coor- 
dinate position of the next point; 

(g) comparing the position difference to an error limit; 

to (h) setting the added distance component (dx. dy, dz! according to a predetermined distance if the position 

difference is within the error limit align he sensor (100) to a new point (23( i napping volume 

along another coc dim , - > ' " 4i i i " m ti i J ^ ^cw point with 

the sensor to establish a new point coordinate position and repeating steps (d) - (g) along the other coordinate 

ts oe component (dx, dy, dz) by decreasing the value of the added distance component 

:: the position difference is not within the error limit and establishing an intern ecii its point (220b) by repeating 
iti ;s (d (g) along the same c rd e axis 

8. The method >f claim 6 im 7. whet ?libt lion method i J 1 re r ping >li ri 

9. The method of claim 8, wherein the mapping volume is approximately 20 cm x 20 cm x 20 cm. 

10. The method of any of claims 8 to 9, wherein the error limit is < 1 mm. 

25 11. The method of any of claims 6 to 10, wherein the sensor (100) is moved by a distance ranging from about 2 cm to 
about 3 cm. 

12. The method of any of claims 6 to 1 1, wherein each time the sensor (100) is moved it is by a constant distance. 

30 13. The method of any of claims 6 to 1 1 , wherein the added distance component is decreased through division by a 
factor of two (X, + dx/2, Y, + dy/2, Z, + dz/2). 

14. The method of any of claims 6 to 13, wherein the sensor (100) is moved according to the vertices (210, 220, 230, 
240, 250) of a cube (200a, 200b). and the entire mapping volume comprises a plurality of cubes. 

35 

15. The method of any of claims 6 to 14, wherein the sensor (100) is moved by a robot arm (1 10). 



Pateritanspriiche 

40 

1. Kaiibrie ! das folgendes umfaRt: 

lehrere Ft - s Feldabstrahier mehrere Abstrahieiernente (32-48) aufweist 

undjede„ 'XL i > i ■> v „ l J j - it> Abstrahieiernente 

45 (32-13) zueinander unterschledlich sind; 

einen Signalprozessor (25), der betriebswirksam mit den Feldabstrahlem (56, 58, 80) verbunden ist; 
©inert PrufpositJonssensor (100) mit dre: Sensorspulen (102, 104, 106); und 
eine Posilionierurigselnrlchtung (110): 

50 fjo i a.iprozessor (25) und m:it der Pos.itionierungs 

einrichtung (110) verbunden ist; 

die -\ ~ * ^ o r " "i betnebswirksam mit r r , g r- verbunden * 

>rpragrammiert ist, 

das ein Abbiidungsvolumen det:nie:1. das :rn wesentllche;-: zu emern Betnebsvolomen ausgenchtet Ist, das von den 
55 Abstrahiern (56 58, 60; erzeugf: wird. wenn diese angeregt weiden; 

der Slg ilpr^ze 0 i tet Ist Sen Prufposiiionssensor (100) zu jeder der Prufpositlonen mittels 

der Positionieru i ig ' 10} zu fuhren; und 

wobei der S a st. ein Abtastsignai von dem Sensor (1 00) zu empfangen, welches 
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das an dem Sensor abge ; < i das Abta lagnetisches 

Feid an dem Sensor definiert, und der Signalprozessor auRerdemdafureingerichtet 1st, eine bekannte theoretische 
magnetische Feldstarke fur jede Prufposition rnit dem tatsachiioh gemessenen Feid zu korrelieren und die Korre- 
- magnetischen Felder zu speiehern. 

System nach Anspruch - be: dem die Posit:onierungseinriohtung elnen Roboterarm (110) autweist. 

System nach Anspruch 1 oder2, bei dem die Sensorspulen (102, 104, 106 c - 1 ueinandersind. 

^ny r ahi u i 3 or bei einei srn (90) nac u c i c y einr <e 

Feid zum Verfolgen einer Position einer Einrlchtung zu erzeugen, wobel das Verfah hi tte aufweist: 

Bereitstellen mehrer Felderzeuger (56, 58, 60), wobei jeoe i - iei aufweist und 

r Al 13' m ' t -rt " ■" 'n\'i t cJ , a-uj^c xmAnregen 

der Abstrahlelemente zueinander unterschiedileh sind; 

Bersfeteiien eines Signalprozessors (25), derbetriebsvvirksam mit den Feldabstrahlem (56, 58, 60) verbunden 
ist; 

Bereitstellen eines Prufpositicnssensor (100) mit are: zueinander orthogonalen Sensorspulen (102. 104. 106): 
und 

Bereitstellen einer Positionierungseinrichtung (1 10), wobel der Prufpositionssensor (100) betriebswirksam mit 
dem Signalprozessor (25) und mit der Positionierungseinrichtung (110) verbunden Ist und die Positionierungs- 
einrichtung (1 10) betriebswirksam mit dem Signalprozessor (25) verbunden ist; 

Vori ^ it is- a .15. i it i- - i f ^ - net ernes Koc-rdsnatencystems, 

das ein Abbildungsvolumen definiert. das im wesenriicaea za dera Betnebsvolumen ausgerichtet ist, das von den 
Abstrahlem (56, 58, 60) erzeugt wird, wenn diese angeregt werden; 
Einrichten des Pri s isors (100) auf eine spezifizierte Prufposition (X 0> Y 0 , Z 0 ); 

Aufnehmen einer Messung des magnetischen Felds der Z-Komponente unEer Verwendung der Z-Spule (1 06) des 
Prufpositionssensors (100); 

j h 11 w >r "r[ il rn r > c 1 ; un , n 11' < 1 - 1 1 r "( u i -*i 1 , " ,\ ) li n>>C(i ; ot> 

(100) an einen Punkt gelangt der zuvor von derZ-Spule (106) eingenommen wurde; 

Aufnehmen einer Messung des magnetischen Felds der Y-Komponente unter Verwendung der Y-Spule (104) des 
Priifpositionssensors (100); 

t h ttwt - Setreiben Posit wg nri< un I 10 o daB die X puis i id jfpositionsse sons 
(100) an den gieichen Punkt gelangt, der vomer durch dis I C ng und durch 

die Y-Spule (104) wahrend der zweiten Messung eingenommen wurde: 

' jfn h r n < r i t I i > \ j n x I " < 

PrOlpositionssensors (1 00); 

Korrelieren bekannter theoretischer Magnetfsldstarken mit den tatsachlich gemessenen Feidern; und 
Spelchern der Korreiationsergebnisse zum Korrigieren der theoretischen magnetischen Feider. 

Verfahren nach Anspruch 4, das weiterhin die Schritte zum Bewegen des Prufpositionssensors (100) zu einer 
nachsten vorbestimmten Prufposition aufweist und VViederhoien o'er Schritte des Aufnehmens der Messungen und 
des schrffiwessen Betreibens der Positionierungseinrichtung (110). 

Kalibrierungsverfahren zum Gebrauch bei einem S stem 90 nach Anspruch 1 das in der Lags t y <e 

Feid zum Verfolgen einer Position einer Einrichtungzu erzeugen, wobel das Verfahren - it> e; etallische Objekte 
berucksichtigi: und foigende Schritte aufweist: 

la) Definieren e:ne ,' Olaidunasaoiumenr i.'JOOa' :nne:halb a'ea err-sugten magnetischen Feicis. 

(b) Plazieren eines metallischen Objekts innema b dt gsvo urnens; 

(c) Ausnchten eines Sensors i 00} an einem ersten Punkt (110) innerhaib des Abbiidungsvoiumens und Messen 
des magneiischen Felds an dem ersten PunKt mix uem Sent -T . t im t rto <c 1 i n * Y 7 
festzulegen; 

id) Bewegen des Sensors an einen nachsten Punkt (220; X; + dx, Yj+dy, Z + dz) ent a iinatenachse 
(210, 220. 230. 240, 250) durch eine hinzugefugte Abstandskomponente (dx, dy, dz) und Messen des magne- 
tischen Felds an dem nachsten Punkt, urn eine nachste Koordinatenposition festzulegen; 
(e) Interpoiieren des magnetischen Felds an einem Zwische' p. - > schen der ersten Position (210) 
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i k n h < v jt i 1 < it (220), ii erne k p l i "a i ■> ' . t n t ulegen. 

(f) Bestimmen des F n J 9n der interpolierten Zwis i tenposition (220a) und 

der tatsachlichen Zwlschenkocrdinatenpositlon; 

tes mit einer Fehlergrenze; 

(h) Setzen (X, Y, Z,) des nachsten Punkts als (Xj = Xj + dx, Y, = Y, + dy, Z = Z; + dz), falls der Positionsabstand 
innerhaib der Fehlergrenze iiegt. und ie i-i der Sch Mite (d; - ;o) ~ hn i eirierarideren Koordiiiatenachse 
{210, 220, 230,0 240, 250}; und 

tzs V ipe -t i t d> u dz) ju - s s c . >rt i 

Abstandskomponente, falls der Positionsabstand nkht innerhaib der Fehlergrenze Iiegt, und Wiederhoien der 
Schritte (d) ■■ {g) entlang der gleichen Kocrdinalenachse. 

7. K 5 rung rfahi « G i ;c beieinei ~ " < ■ i Anspruoh 1 . das in der Lage 1st, ein magro 

Feld zum Verio t on einer Einrichtung zu erzeugen, wobei das Vertahren statische metallische Objekte 

bei-uoksishltg: und foigende Sohrilte aufweist 

(a) Definieren eines Abbildungsvolumens (200a t innerhaib des erzeugten magnetjschen Felds; 
b) Piazieren einer m<-1 i 

(c) Ausrichten eines Sensors (1 00) an einem ersten Punkt (110) innerhaib des Abbildungsvolumens und Messen 
des magnetischen Felds an dem ersten Punkt mlt dem Sensor, run erne erste Koordinatenpositicn (X,. Y.. Z,} 
festzulegen; 

(d) Extrapolieren des magnetischen Felds eines nachsten Punkts (220; X, + dx, Y; + dy, Z; + dz) entlang einer 
Kocrdinateriachse (210. 220. 230. 240. 220; durch erne hinzugefugte Abstandskcmporiente (dx. dy. dz), 

(e) BerechnenderKoordinatenf m i - 3 i j 1 i i ~ i i 
magnetischen reides. urn e:ne extrapokeite Kcordlnatenposltion festzulegen: 

(f) Bestimmen des Positionsabstandes zwischen der extrapolierten Koordinatenpositicn (220) und der tatsach- 
lichen Koordinatenpositicn des nachsten Punkts; 

(g) Vergleichen des Positionsabstandes mit einer Fehlergrenze; 

(h) Setzen der hinzu j« l £ komponente (dx, dy, dz) gema(3 einem vorbestimmten Abstand, falls 
der Positionsabstand innerhaib der Fehlergrenze Iiegt, und Ausrichten des Sensors (100) auf einen neuen 
F 1 * i inn i "i I - n ■■ n v n 1 V 230,240, 
250) und Messen des magnetischen Feids an dem neuen Punkt mit dem Sensor, urn erne Koordinatenposition 
des neuen Punktes festzulegen, und Wiederhoien der Schritte (d) - (g) entlang dci - , u i v 
und 

(1) Setzen der hinzugef A 1 I te (dx, dy, dz) durch Verringem des Wertes der hinzugefugter 

Abstandskomponente falls ue L ^ t ai v uiJ 1 c v all ^ f ^k ui c eg- und Festiegen eines 
Zwrschenpunktes ..:20b: durch Wiedemeien des Schntte id; <g; entlang der gleichen Koordinatenachse. 

8. Vertahren nach Anspruch 6 oder 7, bei dem das Kalibilerungsverfahren fur das gesamte Abbiidungsvolumen ver- 
vollstandigt vvird. 

S. Verfahre 

10. Vertahren nach einem der AnsprGche 0 bis 9. be: dem d:e Fehlergrenze <1 mm betragi 

11. Vertahren nach einem der AnsprGche 6 Ok 10. bei dem der Sensor ( 200; urn einen Abstand bewegt wsrd, der j m 
Bereich zwischen etwa 2 cm und ekva 3 cm iiegt. 

1 i. i ihi t i nil -i 1 1 i 1 i 1 i Sc > ird. dies ui 

einen konstanten Abstand erfoigi. 

13. Vertahren nach einem der Anspruche 6 bis 1 1 J rte durch Division ui 
einen Faktor Zwei (Xj + dx/2, Yj + dy/2, Z, + dz/2) verringert wird. 

14. Vertahren nach einem der Anspruche 6 bis 13, bei dem der Sensor (100) gemaB den Eckpunkten (210, 220, 230, 
240 250) ein <ubt C 1 I bewegt wird und das ges ite At: Cuben ut 

15. Vertahren nach « v ie 6 bis 14, bei dem der Sensor (100) von einem Pv be 
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Reversdications 

1. Systerne de calibrage (90) comportant : 

■ une pluralite de generators de champs {56, 58, 60), chaque generateur de champ ayant une pluralite 
d'eiements :omiantaerieiateursi32a48). cheque eiemenf:omiai:tgenerateu:-p;ode:santun champ magnetique 
distinct d'un autre par une mise sous tension s mi. I 32 a 48), 

■ un prooesseurde signal (25 e e ae r ■> n - v i 3 ^3 60) 
K un capteur de position de test (100) ayanttrois bobines de capteurs (102, 104, 106), et 

m un dispositif de positionnement (110), 

M dans lequel ledit capteur de position de test (100) est relie de maniere fonctionneiie audit processeur de 
signal (25) et ii ie f onnement (1 10), 

^ dans leque! ledit dispositif de positionnement (110) est relie de maniere fonctionneiie audit processeur 
de signal (25), 

■ dans lequel ledit prooesseurde signal (25) est pre-programme en utihsant une pluraiite de positions de 
test dans ti i tnn>. l .ant i< < d I > i >o ^ 
volume de fonctionnement produit par Iss generateurs (58, 58. 60) lorsque ceux-ci sont mis sous tension, 

■ dans lequel ledit prooesseurde signal (25) estadapte pour guider ledit capteur de position detest (100) 
jusqu'a chacune desdites positions de test au moyen dudit c p - tionnement (1 10), et 

M dans iequel iedit processeur 3 1 j it i i ic provenant 

audit capteur (100), Indicatif du champ lagneti t n I o fe detection 

deflnissant un champ magnetique mesure sur ledit capteur, ledit processeur de signal etant egalement 
adapte pour correler une valeurd ie pourc \ etestave 

le champ mesure reel, et pour memoriser des resultais de correlation pour une correction des champs 
magnetiques theoriques. 

2. Systerne selon la revendication 1, dans iequel ie dispositif de positionnement comporte un bras de robot (1 10). 

3. " i i i i*eurs (102, 104, 106) sont mutuellernent 
orthogonaies, 

4. Precede de caiibrage pour une utilisation avec un systerne (90) selon la revendication 1 capable de produire un 
champ magnetique pour suivre une position d'un dispositif, le procede comportant ies etapes consistant a : 

■ agenceruneplur ieciege i& champs (56. 58. 60}, chs a ant nents formant 

i j > (32 a 48), chaque element formant generateur ef anl adapte pour produire un champ magnetique 
t i-pi ^ i - u n I - cnt„ i -ntt it- 

■ agencer un processeur de signal (25) relie de maniere fonctionneiie auxdits generateurs de champs (56, 58, 
60), 

■ agencer un capteur de position de test (100) ayant trois bobines de capteurs mutuellernent orthogonaies 
(102, 104, 106), et 

& agencer un ciispooitit de positionnement ( '.'■(;). dans lequel iedit capteur cie position de tost { 1 00) est relie do 
nan ere to iior is j audit pre e .- s u v . ~ . . t 1 . . ' er "0 t c 2 lis x -it 

de positionnement ( 1 10) est relie de maniere fonctionneiie audit processeur de signal (25), 
m pre-prog irmnpr edi processeu de ig mtun i i u ten 

3 nnee ieirni s neoec i<= ^ti ec un volume de fonctionnement 

produit par Ies generateurs (56, 58, 60) lorsque ceux-ci sont mis sous tension, 
U amener le capteur ce t v ti n v c "0 , ^ _ i Jc 1 | ' , , 

■ prendre une mesure du champ magnetique de la composanie Z on utiii; anl i. • bobine Z. I Of. • du capteur de 
position de test (100), 

■ faire progresser pas a pas le dispositif de positionnement (110) de sorie que ia bobine Y (104) du capteur 
v x ,i c 1 de *e~t J .0 r < , ^ ■> . c emment occupe par ia bobine 7. (106), 

■ prendre une mesure du champ magnetique de la composanie Y en utiiisant la bobine Y (104) du capteur de 
position de test (100), 

■ faire progresser pas a pas le dispositif de positionnement (110) de sorte que la bobine X (102) du capteur 
de position de tost (100) est sur ie meme point precedemment occupe par la ore Z (106) pendant la premiere 
mesure et par ia bobine Y (104) pendant ia seconde mesure, 
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m prendre une mesure du champ magnetique de lacomposante X en utiiisant la bobine X (102) du capteur de 
position de test (100), 

^ correler des valeurs de champs magnetiques theoriques connues avec les champs rnesures reels, et 
^ memorise!' des ies^ t i J - i i c i pour une correction des champs magnetiques theoriques. 

Precede selon ia revendication A, comportant egaiement les etapes consistant a den de position de 

test (100) jusqu'a une position de test predetermines suivante, et repeter ies etapes consistant a prendre des 
mesures et a faire progresses pas a pas le dispositif de positionnement (110). 

Precede de calibrage pour une utilisation avec un systems (90) selon la revendication 1 capable de produire un 
champ magnetique poursuivre une position d'un dispositif, le precede representantdt t > % 

> 1 comportant les > 1 ,a consistent a : 

(a) definir un volume de cartographie 2 ' >r :> magnetique produit, 

( ^ er 1 c de \ to s \y e 

(c) aligner un capteur (100) sur un premier point (210) a ilnteneur du volume de cartographie et mesurer le 
champ rnagaeiique sur le premier point a I'aide du capteur pour etablir une premiere position dans ie systems 
de coordonnees (X M Yi, Z;) , 

id) depiacer le capteur jusqu'a un point suivant (220 : X; + dx, Y, + dy, 2, + dz) ie long d'un axe de coordonnees 
(210. 220, 230 240. 250) d une vaisur ©gale a une ccmposante de distance ajoutee (dx, dy, dz) et mesurer ie 
champ magnetique sur !e point suivant pour etablir une position suivante dans ie systeme de coordonnees 

(e) interpoier ie champ magnetique sur un point intermediate (220a) entre la premiere position (210) et la 
position suivante dans le systems de coordonnees i c x c - x une position intermediaire interpolee dans 
ie systeme de lee 

(f) determiner ia difference de position entre la position i it xmedia i e inter poiee dans ie systeme de coordonnees 
(220a) et la position intermediate reelle dans le systeme de coordonnees, 

(g) comparer ia difference de position a une llmite d'erreur, 

(h) fixer (X,, V,. Z,> du point suivant a (X, - X, + dx, Y, - Y, + dy, Z, - Z, + dz) si ia difference de position est dans 
ia limite d'erreur et repeter les etapes (d) a (g) le long d'un autre axe de coordonnees (210, 220, 230, 240, 250), et 

< a a r c x i u i u 1 ia a 3 t J 

ajoutee si la difference de position n'est pas dans la llmite d'erreur et repeter ies etapes (d) a (g) ie long du 
meme axe de coordonnees. 

Precede de oaiibrage pour une utilisation avec un systeme (90) selon la revendication 1 capable de produire un 
champ magnetique poursulvre une posifion J disposifif. ie precede repreaerifanf des objets rnetaiiiques sfaiiques 
et comportant les etapes consistant a ; 

(a) definir un volume de cartographie (200a) dans le champ magnetique produit, 
c . , x \ ' ' , "I 1 , ->> de cartographie, 

(c) aligner un capteur (100) sur un premier point (210) a i'mterieur du volume de cartographie et mesurer le 
champ magnetique sur le premier point a I'aide du capteur pour etablir une premiere position dans le systeme 
de coordonnees (X,, Yi, Zj), 

(d) extrapoler le champ magnetique d'un point suivant (220 : X, + dx, Y, + dy. Z, + dz) le long d'un axe de 
joorcfonnees (2 20. 230. 240, 250) d eui ie distance sjoutee (dx dy dz) 

(e) calculer la position dans ie systeme de coordonnees sur le point suivant extrapoie (220) sur ia base du 
champ magnetique extrapoie pour etablir une position exirapoiee dans ie systeme de coordonnees, 

(f) deten iner |[ e pc ition e * e a exirapoiee c - a j e c nees (220) et 
la position reelle dans ie systeme de coordonnees du point suivant, 

(g) comparer ia difference de position a une simile d'erreur, 

dz) on fonction d'une nines si la difference 

de position e 1 dane a lln z d'erreui aligr e ; capteu 00) sur un nc c 1 30) a iterieur du 5 u i c 

de cartographie ie long d'un autre axede coordonnees (21 0, 22f 30 40, 250) et me ; ha p magnet i 
sur ie nouveau point a I'aide du capteur pour etablir une position du nouveau point dans le systeme de coor- 
donnees et repef - (d) a (g) e ong de I'autre axe de coordonnees, et 

(!) fixer ia composante de distance ajoutee (dx, dy, dz) en reduisant la vaieur de la coniposanfe de distance 
ajoutee si la difference de position n'est pas dans la limite d'erreur et etablir un pc nl interrm t&aire :220b > en 
repeiant les etapes (d) a (g) ie long du meme axe de coordonnees. 
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8. Precede selon la revendication 6 cu 7, dans iequel le precede de calibrage est complete pour la totalite du volume 
de cartographie. 

9. Proced<- d ration 8 dans Iequel le volume de carte .i jph , i - ' r ent de 20 cm x 20 cm 

x 20 cm. 

10. Precede selon I'urte quelconque des revendications 6 a 9, dans Iequel la limite d'erreur est < 1 mm. 

11. Procede selon I'une quelconque des revendications 6 a 1 0, dans Iequel ie c - 
allant d'environ 2 cm jusqu'a environ 3 cm. 

12. Procede selon I'une quelconque des revendications 6 a 11. dans Iequel a chaque fois que le capteur (100} est 
depiace, ii Test d'une distance constants. 

13. Procede selon I'une quelconque des revendications 6 a 11. dans jequei la composante de distance a:outee eat 
reduite par une division par yr> raoteur de deux (X, + dx/2, Y, + dy/2, Z, + dz/2). 

14. Procede selon I'une quelconque des revendications 6 a 13, dans iequel le capteur (100; est depiace en fonction 
J' r " " 0 ° J b 1k 5 U! i l i 

une pluralite de cubes. 

15. Procede seion I'une quelconque des revendications 6 a 14, dans Iequel ie capteur (100) est depiace par un bras 
de robot (110). 
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